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Abstract 


The resin transfer molding (RTM) manufacturing process consists of either of 
two considerations. The first is the fluid flow analysis where tracking of the flow 
fronts are of primary importance, and the second is the combined multi¬ 
disciplinary flow/thermal analysis, where kinetics and cure need to be 
accounted for in the process modeling stages. In the combined multi¬ 
disciplinary flow/thermal analysis, viscosity is a function of temperature. The 
so-called continuous sensitivity formulations are developed for the latter case of 
non-isothermal process modeling of composites manufactured by RTM to 
predict, analyze, and optimize the manufacturing process. Attention is focused 
on developments for non-isothermal process modeling simulations, and 
numerical illustrative examples are presented for the sensitivity analysis of 
structural composites manufactured by RTM. 
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1. Introduction 


Process modeling employing resin transfer molding (RTM) to manufacture complex struc¬ 
tural geometries involves the injection of a polymer resin into a porous fibrous preform 
conforming to the geometry of the manufactured part. The mold is normally at a higher 
temperature than the resin; therefore, the analysis of this process requires accounting for 
both the flow and thermal-kinetic equations. The solutions of these equations provide the 
flow front position, temperature, and cure information to the analyst. After the results from 
the numerical analysis are obtained, a systematic process of optimization is performed with 
the assistance of the continuous sensitivity equation, thereby avoiding heuristic trial-and- 
error methods. 

The thermal processes involved in the resin transfer molding (RTM) model are investi¬ 
gated in this report. The non-isothermal RTM analyses have been previously investigated in 
references [1-4]. Modeling of both temperature and cure are important because the viscosity 
of the resin is sensitive to both of these factors. Before reaction of the polymer resin occurs, 
the resin is temperature sensitive; but after polymerization begins, the resin is cure sensitive 
[1]. Although the continuous sensitivity equation (CSE) has been previously applied under 
isothermal conditions [5,6], the aim of this report is to extend this to non-isothermal situ¬ 
ations for the first time. Besides the development of the CSE for non-isothermal RTM, the 
fill time sensitivity with respect to the inlet temperature is also described. 

The CSE for mold filling is formulated for process modeling of composites that take 
into account non-isothermal considerations. The thermal model is first presented and solved 
numerically by coupling with the flow model. This thermal model includes solutions for 
temperature and cure. The CSE for non-isothermal RTM filling is presented next, followed 
by the CSE for the temperature and cure sensitivity considerations. Finally, illustrative 
numerical results are presented for two sample axi-symmetric geometries, which includes 


1 



verification of the fill time sensitivity for a sensitivity parameter of inlet temperature. 


2. Non-Isothermal Resin Transfer Molding 


In this section, the temperature and cure model equations are described for non-isothermal 
RTM filling. Non-isothermal effects are an important consideration in many manufacturing 
process modeling applications employing RTM. Heat is transferred via various mechanisms 
of transport during a mold filling process. These include convection from the mold walls, ad- 
vection from the fluid, radiation to and from the mold, and heat generated by the exothermic 
curing process of the resin. For non-isothermal considerations, the independent parameters 
include conductivity, permeability, viscosity, and inlet pressure or flow rate. The dependent 
variables are the temperature and amount of resin cure. 

The first step to developing the so-called continuous sensitivity equations for the non- 
isothermal filling process is to begin from the governing model differential equations. By 
assuming that the fiber and resin temperatures are at the same equilibrium temperature 
inside a small control volume [1], the energy balance can be written as 

pc p ^ + Pr c pT (u • VT) = V • kVT + §G (1) 

where G is the heat generated from the curing of the resin and $ is the volume of fiber 
fraction. The average material properties are computed by 

pCp — Q PfCpr ~\~ (1 ^^PfCpf 

k = k s + k D (2) 

k s = $k r + (1 — &^kf 

where the subscripts / and r denote fiber and resin properties, respectively. k D is the 
thermal dispersion conductivity due to mechanical dissipation. The computation for the 
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effective conductivity is given as 



kfk r 

kfW r + k r Wf 


(3) 


where 

— w r 
JL 

_ Pf 
' r ± , (l-») 

Pf Pr 


Wf = 1 


(4) 


Equation (1) represents the equilibrium temperature model used for developing the finite ele¬ 
ment equations that can then be readily implemented and solved numerically. The boundary 
conditions for the non-isothermal mold-filling process are given as 

T = T w at mold wall 
T = T r0 during filling, at inlet 

dT ( 5 ) 

k — = (1 — §)p r c pr u • n(7/o — T) at resin front 

T(t = 0 ) = T W 


where T r is the temperature of the resin, T w is the wall temperature, and n is the direction 
normal to the resin flow front. The initial conditions-for solving the temperature problem 
are given as 

T r (t = 0) = T r0 and Tf(t = 0 ) = T W (6) 

where Tf is the temperature of the fibrous preform. Another thermal computation involved 
in non-isothermal RTM filling is cure. The curing analysis begins with the species mass 
balance and is given as 

r\ 

+ u ..Va = $Jk (7) 
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Here, a is the degree of cure, u is the velocity field, and R a is the rate of chemical reaction. 
The boundary conditions for the cure problem are given by 

a = 0 during filling, at the model inlet(s) (8) 


The initial conditions are given by 

a(t = 0) = 0 (9) 

Since curing is an exothermic process, the heat generated is calculated with the following 
equation [2-4] 

G = H R R a (10) 


where the term R a is defined as 

R a = (K x + K 2 a nl )(l - a) n2 


( 11 ) 


and Hr is the heat of reaction per unit volume for the pure resin. The constants K x and K 2 
are assumed as 


K\ = A x exp 


K 2 = A 2 exp 



( 12 ) 


where A x , A 2 , E%, E 2 , n x , and n 2 are the kinetic constants determined experimentally for 
each resin used. The resin viscosity is a function of the degree of cure and temperature. The 
model used for the viscosity is given as [7] 

" = A “ exp (It) (5^) (13) 

which is used in conjunction with the implicit filling technique due to Mohan et al. [8]. This 
method has proven to be an effective approach for solving this class of problems. 
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Continuing with the development of the heat transfer finite element equations, begin with 
the thermal equilibrium model in equation (1) and employ the method of weighted residuals. 
Thus, 

dQ = 0 (14) 


L 


w 1 


' dT 

pc p — + p r c pr (u • T) - V • kWT - 


Separating the terms and integrating by parts yields 

[ W T pc p ^~ dQ + [ W T p rCpr {u ■ VT) dtt - [ W T (kWT ■ n) dT+ 
Jvt ut J n J r 

I 


VW T kVTdCl- [ W T $GdQ = 0 


(15) 


Defining the weighting functions W to be the Streamline Upwind Petrov-Galerkin (SUPG) 
weighting functions N + A(u • VN) and interpolating for T yields 

W = N + A(u • VN) 

num. nodes ( 16 ) 

T = 

1=1 


where i represents the associated node numbers. The heat flux, q, is defined as 

fcVT • n = —q • n 


(17) 


Substituting equation (16) and equation (17) into equation (15) yields 

[ Wipc p Ni dtt] ^+\ f W lPr c pr (u • VNi) dQ Ti+ [ VWikVNidQ T< 

.Ja J ot [J a J Un J ^g\ 

= [ Wi(kVT-n)dV+ [ Wi$GdQ 

Jr Jn 


For the temperature analysis, the resulting semi-discretized equation is thus obtained as 


CT + (K ad + K cond )T = Q q + Q 6 


(19) 


o 




where the terms in equation (19) are defined as 

C = [ W T pCpN<Kl 
J n 

K ad = [ WVvCu-B/v) dQ 

Jn 

Kcond = [ dQ • 

Jo. 

Q q = J^w T (-q-n) dr 

Q 6 = [ 

7n 

The time discretization for the temperature analysis is employed, for an arbitrary 9 , as 


T„ = 


“ 71+1 -*-71 


The cure model equation is also solved using the finite element method. Applying the method 
of weighted residuals on the species mass balance in equation (7) yields 


J W T + u -Vq- $R a J dQ = 0 (22) 

Subsequently, separating terms and integrating by parts yields 

f dQ + f W r u • VadQ - [ $R a d<R = 0 (23) 

Jn ot Jn Jn 

As in the temperature problem, the weighting functions W are defined as the SUPG weight¬ 
ing functions N + A(u • VN) and interpolating for a yields 

W = N + A(u • VN) 

num. nodes (24) 

a= WiCti 


Rewriting equation (23) in the semi-discretized form results in 

Cd + Ka = Qr q 
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where the terms in equation (25) are defined as 

C = [ $W t N dQ 

Jn 


K 


= [ W T (u • Bjv) 
Jn 

Qr q = f $w T R a dn 

Jn 


dQ, 


(26) 


The time discretization for the cure analysis is employed, for an arbitrary 6, as 

&n +1 


OLB = 


At 


(27) 


A more detailed discussion of the non-isothermal RTM computational analysis can be found 
in Ngo et al. [9]. 


3. Sensitivity Analysis of Non-Isothermal Resin Trans¬ 
fer Molding 


In this section, the continuous sensitivity equation for non-isothermal RTM filling, tempera¬ 
ture, and cure are defined. The computational procedure for the non-isothermal RTM filling 
problem is described, along with some illustrative results from an example axi-symmetric 
model. Finally, the fill time sensitivity with respect to the inlet temperature is discussed 
and verified with the use of an additional axi-symmetric example. 


3.1 RTM Filling Sensitivity 

For the non-isothermal pressure sensitivity calculation, viscosity is a function of temperature 
and cure (i.e., n (T, a) ). The sensitivity equation for RTM filling is given by the following 
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semi-discretized equation as 


-w—P + KS P = S q 
dp 


(28) 


where p is the specified sensitivity parameter. Equation (28) is the same CSE used in isother¬ 
mal filling simulations and the derivation thereof is given in the aforementioned isothermal 
references [5,6]. The terms ^,K, and S q from equation (28) are defined as 



In equation (29), K is the permeability tensor of the fiber preform. With the viscosity defined 
as a function of temperature and cure, this leads to an additional term in the non-isothermal 
RTM filling continuous sensitivity equation, namely ^ (jffiaj') 

£* = [ B T^Bd(l+ f b t k (- 4 ) ^ ( J’ a) B dtl (30) 

dp J dp p J V M 2 / dp 

n n 


where d ^ ~ is defined as 

dfi(T(p),a(p);p) _dp . 
dp dp * 


a„ 


a g -a 


A+B a 


E„ 


,(fr) 


+ A^e 


(It) 


CXn 


CZq 


a 


RT 2 

A+B a r A + B a 


Ota 


a 


dT 

dp 

+ Bln 


Ota 


a g — a 


da 

dp 


(31) 


Note, that for p = constant and §£ = 0 where p ^ p, it will yield the same results as the 
isothermal model equation. In the non-isothermal filling simulations the temperatures are 
also computed, which allows for the solution of the temperature sensitivity equations. 
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3.2 RTM Temperature Sensitivity 


The temperature sensitivity is evaluated by taking the partial derivative of equation (1) with 
respect to the sensitivity parameter p. 


d_ 

dp 


dT 

pC p — + p r Cp T {U • VT) 


d_ 

dp 


V-A;VT + 




This results in 


| w £ ++ £<*><“' VT)+ ^| (u ' vr) 

= V • ^ VT + V • kVS t + 

ap ap op 


(33) 


where the temperature sensitivity is defined as St — Moving all the terms in equa¬ 
tion (33) to the left hand side and applying the method of weighted residuals yields 



!<*»£ + + | to) (■ ' VT) + (u . VT) 


„ dk„„ „ 

V • -VT + V • A:VS r + w-G + 
ap ap ap 


dfi = 0 


dQ— 


(34) 


Before deriving the finite element equations, the partial derivative of the boundary con¬ 
ditions for the temperature analysis need to be computed. For the temperature analysis, 
equation (33), the boundary conditions are as follows: 


St 


St — 


dT w 

dp 

dT r0 

dp 


at the mold wall 


during filling, at mold inlet 


(35) 


k'S^- = §)p T c pT (u • n)(T/ 0 - T)) at resin flow front 

on dp 


After applying the Green-Gauss theorem to convert the volume integral to a surface integral, 
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the temperature sensitivity equation is given as 


L wT ! to) f dn+ i + L wT r P tov) («• * 

+ J WV*> f0.r) <KJ + jf WV (u • VS r ) dfi - jf W T (|jvr • n) dr 

+ f VW T • — VTdQ - [ VW T (kVS t -n)dT+ [ VW T • W5 r dfi 
Jn dp Jr J n 

- f W T ^-GcM - [ W r $^ dfi = 0 
Jn dp Ju dp 


where pc p is defined in equation (2) and p r c pr are the resin material properties. Defining 
that the weighting functions W be the SUPG weighting functions N + A(u • VN), and 
interpolating for T and Sj yields 

W = N + A(u • VN) 

num. nodes 

T = E W ‘ T ‘ (37) 

2=1 

num. nodes 

St = E W ‘ S n 

2=1 


Substituting equation (37) into equation (36), yields the semi-discretized finite element equa¬ 
tions for the temperature CSE, given as 

~T + CS T + (^ T + (K ad + K„„d)S T = S q + S 6 (38) 

dp \ dp dp J 
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where each of the terms in equation (38) are defined as 

f = /. 

C-[ W T (pep) N dQ 

Jq 

= i wT i (u - B)dn+ L wr ^ (S' B ) dQ 


dK 


K ad = [ W T (p r C pr ) (u • B) dQ 
Jn 

dp ~ So. 




dfce// 

<9p 


B dQ 


K C ond — / B|y/^gyyBdQ 

do 

s q = jf w r VT • n) dr + ^W T (k\/S T ■ n) dr 


(39) 


with G defined in equation (10). The time discretization is employed as 

Tn+i T n 


T„ = 


St 9 — 


At 

^Tn+l ^Xn 

At 


(40) 


The fully discretized representations are thus given as 


[C + 0KAt] S Tn+1 = — 


dC a dK A ' 

La? + s a? Ai 


■ n+l 


+ 




[dp 


dp 


+ [C - (1 - 0)KAt] S Tn + At 


(l-«)^+/ F " +1 


dp 


dp J 


(41) 


with temperature T and temperature sensitivity St approximated by 

T* = (1 — d)T n + dT n+1 
St # = (1 - ^) S T„ + ^S Tn+1 


(42) 
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and F is defined as 


(43) 


F = S„ + S, 


Some of the terms necessary for calculating the temperature sensitivity need to be defined. 
The first comes from equation (3), where is defined from the following equation for 
sensitivity parameter p = k r , 


dk 


eff 


ktk 


fro r 


dp kfW r + k r Wf (kfW r + k r Wf) 


;Wf 


(44) 


Or for p = kf, 


dk, 


'eff _ 


kfk 


fli r 


dp kfW r + k r Wf ( kfW r + k r Wf) 


rW r 


(45) 


The heat generation term, is defined as 


dG _ „ dR a 
dp Hr dp 


(46) 


where Hr is the heat of reaction per unit volume for the pure resin. The rate of reaction 
term, is defined by 

- rv '*" 2 

dp' 

,da 


dp a { dp dp 


a" 1 - + K 2 n 1 a ni -r-)(l - a) r 


(47) 


+ (K 1 +K 2 a m )n 2 (1-q) ( ” ! - 1 > 


dp 


where and ^ are found by taking the partial derivative of K\ and K 2 in equation (12) 
with respect to the sensitivity parameter p, i.e., 


dK x E X A i ,_^)dT 

- — -gv rt ) - 


dp RT 2 


dK 2 E 2 A 2 iiz)dT 

- — - RT J 


dp RT 2 


dp 

dT 

dp 


(48) 


Equation (41) is used to solve for the temperature sensitivity profile for each time step after 
the temperature results have been calculated. 
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3.3 RTM Cure Sensitivity 


In the non-isothermal RTM analysis, cure is also computed. This leads to the following 
evaluation of the cure sensitivity equations. The cure sensitivity equation is computed by 
taking the partial derivative of equation (7) with respect to the sensitivity parameter p, i.e., 

£(*| + u -v<»=**.) (49) 


This leads to 


d$da S a dn _d$ 

dp dt dt dp dp 


+ $ 


dR a 

dp 


(50) 


where the cure sensitivity is defined as S a = Applying the method of weighted residuals, 
yields 


/ 


W 1 



+ ^ + %'^ a + u ' VS ° d P Ra $ dp 


dQ, = 0 


(51) 


As in the temperature problem, the weighting functions W are defined to be the SUPG 
weighting functions N + A(u • VN), and interpolating for a and S Q yields 

W = N + A(u*VN) 

num. nodes 

a = Wi0ii (52) 

2=1 

num. nodes 

s„ = y , w < s « 

2=1 


The semi-discretized cure sensitivity equation is thus given as 


dC . <0K ^ _ 5Qr q 

—a + CS„ + -g^a + KS„ - dp 


(53) 
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where the terms in equation (53) are defined as 


C 

dC 

dp 

K 

dK 

dp 

3Qr q 

dp 


f $W r N dtt 

[fw-NiO 

Jn dp 

[ W t (u-Ba r )dtt 
Jn 

,du 


= [ W t (^-B N )dQ 
Jn dp 

= f T R a dtt+ [ $W- 
Jn dp Jn 


- dR a 
dp 


d£l 


with the time discretization employed as 


CX0 = 



^n+l ^-n 


At 

S Qn+1 — S Qn 
At 


The cure sensitivity is solved numerically employing 



IdC n dK ' 


dC „ 

dK 1 

[C + OKAt] S Qn+1 = - 

"X—1“ 

dp dp 

CX-n +1 + 


- e) -§i ht \ 


+ [C - (1 - 0)KAt] S Qn + At 




+ e 


dF~ n+l 

dp 


with cure a and cure sensitivity S a approximated by 


olq = (1 — 0)a n + 9oc n+ 1 

Sag = (1 — ^)S Qn + 0S Qn+] 


All of the terms in equation (56) are defined in equations (26) and (54). 



3.4 Computational Procedure 


The 14 steps required to perform the sensitivity analysis for non-isothermal RTM process 
modeling are outlined next. 


1. The filling analysis is described by the following semi-discretized equation: 

KP = q (58) 

2. For non-isothermal filling analysis include the necessary temperature, equation (19), 
and the cure equation, equation (25). 

3. The CSE for the RTM filling simulation is given in semi-discretized form in equa¬ 
tion (28). 

4. For the non-isothermal sensitivity analysis, it is also necessary to take the partial 
derivative of the temperature and cure equations with respect to the sensitivity pa¬ 
rameter p, equations (33) and (49). This results in the semi-discretized equations for 
temperature and cure, equations (38) and (53), respectively. 

5. To solve the filling equation, the boundary conditions are defined as 

dP n ., „ 

—— = 0 on mold walls 

an 

P = 0 at flow front 

P = P 0 prescribed pressure at inlet (59) 

or 

q = q 0 prescribed flow rate at inlet 

6. To solve the filling sensitivity equations, the boundary conditions need to be computed 
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by taking the partial derivatives of all filling boundary conditions given as 


— (!¥L'\ = 0 on Mold Walls 


dp \dn 
dP 


= 0 at flow front 
dp 

dP dP 

— = —- for constant pressure at inlet 
dp dp 


or for constant flow at the inlet 


( 60 ) 


— = for constant flow rate at inlet 
dp dp 

7. For the non-isothermal analysis, the temperature and cure boundary conditions are 
required and are defined in equations (5) and (8), respectively. 


8. For temperature and cure sensitivities, the thermal boundary conditions need to be 
computed. The temperature sensitivity boundary conditions can be found in equa¬ 
tion (35), with the cure sensitivity boundary conditions defined at the inlet as S a = 0. 

9. For the non-isothermal RTM sensitivity analysis, ^ and S g must be computed and 
boundary conditions applied. The matrix ^ requires ^ and " as defined in 
equation (31). The values in vector S g are set to zero, except where the inlet flow rate 
has been defined. 


10. The temperature sensitivity analysis requires computation of five matrices and vectors 

including, fg, 9K ^° nd , S q , and S 6 . The values required for all of these terms can 

be found in equation (39). 

11. Similarly, for the cure sensitivity, fS, f|, and are required and are defined in 
equation (54). 

12. Once all of these equations are defined and the boundary conditions applied, the results 
must be obtained. First the pressure results are computed, followed by the pressure 
sensitivity results. After the filling is completed for the current time step, the thermal 
computations are performed. The temperature and cure considerations are obtained 
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first, followed by the temperature and cure sensitivities. Note that the temperature 
and cure time step is normally less than the filling time step (see Ngo et al. [9]). 

13. Before continuing on to the next filling step, compute the new viscosity and viscosity 
sensitivity values from the thermal results and reform K and 4^. 

14. Proceed to the next filling time step until the mold is completely filled. 

3.5 Example Geometry With Numerical Results 

An axi-symmetric example geometry is shown in Figure 1. The disk is filled with inlets 
placed along the inside radius of the disk. By varying the inlet boundary conditions and 
the material properties, it is possible to show how these variables affect the fill time or inlet 
pressure of the model. Observation of the results of the sensitivity values for these variables 
are used to indicate the qualitative and quantitative changes these variables produce. The 
results presented in Figures 2-5 are obtained with the geometry and mesh shown in Figure 1. 
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Fill Time 



(a) 4.00-in. R disk with 1.00-in. R hole for non- (b) Mesh of 4.00-in. R disk with 1.00-in. R 
isothermal RTM study hole 

Figure 1: Example geometry and finite element mesh. 




(a) Non-isothermal fill time vs. inlet pressure (b) Non-isothermal fill time sensitivity vs. in¬ 
let pressure 


Figure 2: Non-isothermal fill time and fill time sensitivity vs. inlet pressure plots for the 
4.00-in. R disk model. 
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(a) Non-isothermal fill time vs. inlet flow rate (b) Non-isothermal fill time sensitivity vs. in- 

let flow rate 

Figure 3: Non-isothermal fill time and fill time sensitivity vs. inlet flow rate plots for the 
4.00-in. R disk model. 



(a) Non-isothermal inlet pressure vs. inlet flow (b) Non-isothermal inlet pressure sensitivity 
rate vs. inlet flow rate 

Figure 4: Non-isothermal inlet pressure and inlet pressure sensitivity vs. inlet flow rate plots 
for the 4.00-in. R disk model. 
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(a) Non-isothermal fill time vs. inlet temper- (b) Non-isothermal fill time sensitivity vs. in- 
ature let temperature 


Figure 5: Non-isothermal fill time and fill time sensitivity vs. inlet temperature plots for the 
4.00-in. R disk model. 


Figure 2 shows plots of the fill time and the fill time sensitivity versus inlet pressure. 
These results demonstrate that the fill time decreases with increasing inlet pressure, which 
follows the trend observed in the analytical solution for isothermal considerations as given 
by 


. 

kP 0 


r R 2 . , „ . - 

~ ln — + 


R\ R 2 R 0 2 


(61) 


where 


H = viscosity 
$ = porosity 
k = permeability 
P 0 = Inlet Pressure 
R 0 = Inner Radius 
R = Outer Radius 
t = Time to fill from Rq to R 
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Figure 3 shows the fill time and the fill time sensitivity versus inlet flow rate. The results 
are as expected since the volume of the mold remains constant; and, the fill time decreases 
as inlet flow rate increases. Figure 4 is a plot of pressure at the mold inlet with respect 
to the inlet flow rate; it shows that inlet pressure is approximately a linear function of 
the inlet flow rate for the given range of values. Figure 5 shows results unattainable from 
isothermal analysis. From the numerical results it is possible to plot the fill time and the 
fill time sensitivity versus inlet temperature. As expected, the fill time decreases as the 
inlet temperature increases because viscosity is a function of temperature, and the viscosity 
decreases with increasing temperature. Hence, from the isothermal analytical solutions, one 
observes that as the viscosity decreases for a constant inlet pressure, so does the mold fill 
time. All of these results follow what is expected from experience and from representative 
isothermal analytical solutions. The sensitivity results also follow the expected trends for 
the specified boundary conditions. 


3.6 Verification of Non-Isothermal RTM Sensitivity Equations 


Since there does not exist an available analytical solution for the non-isothermal filling pro¬ 
cess, the following definition of the derivative is employed: 


,. tfui(p + Ap) — tfui{p) _ dtfui 

inn 7 o 

A p op 



(62) 


The associated error calculation is taken as 

Ap • S tfill 


Error = 


tmip + Ap) -tfm{p) 


- 1.0 


(63) 


where p and Ap are the value of the sensitivity parameter and the change in the sensitivity 
parameter, respectively (e.g., T inlet ). For verification purposes, the error, as defined in equa¬ 
tion (63), should monotonically approach zero as Ap decreases. Also, the results obtained 
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from a Taylor series expansion, equation (64), should closely follow the numerical results for 
small changes in the sensitivity parameter, or graphically these results should be tangent to 
the solution of the numerical analysis obtained by employing repeated numerical simulations. 

The first order Taylor series expansion is given as 

w{p) = w{po) +~\ vo {p-po)+■■• (64) 

where w(p) is the numerical result and p is the sensitivity parameter. Note that the non- 
isothermal model verification in this section refers to the axi-symmetric model shown in 
Figure 6(a), with the following parameters: 

Tmou = 450.00°F 
At = 0.01s 
k = 5.300 • 10 _6 m 2 
Po — 300.00psz 

Since the focus of this work is non-isothermal RTM process modeling sensitivity, the veri¬ 
fication is only performed here for the thermal sensitivity parameter of inlet temperature. 
For the verification of the inlet temperature sensitivity parameter, the Taylor series is used 
to compare the sensitivity results to the derivative of the fill time with respect to the inlet 
temperature. The original Taylor series assumed for estimating the fill time results is shown 
here as 

tfill 1 = tfM 0 + ‘ AT i n i e t (65) 

In the first-order Taylor series estimation from equation (65), the resin viscosity is not 
explicitly included. The viscosity term is instead encapsulated inside the fill time sensitivity 
term, j% fm , in equation (65). Viscosity is a function of the temperature profile inside of 
the part, which in turn is a function of the inlet temperature. This dependence of the 
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(a) Axi-symmetric part (b) Axi-symmetric part mesh for verification study of non- 
for verification study of isothermal RTM 
non-isothermal RTM 

Figure 6: Geometry and finite element mesh used for analytical/numerical comparison. 

temperature profile on the inlet temperature gives the viscosity term a reliance on inlet 
temperature which is used to define a new Taylor series dependent on viscosity. In order 
to explicitly list the viscosity in the Taylor series expansion, the following modified Taylor 
series is used, which is based on the viscosity sensitivity. 

tfuh = tfuio + ^ ( 66 ) 

The A// term is computed for every filled node in the numerical model at each filling time 
step. The final A/x used in the Taylor series is a weighted average of each time step, weighted 
by the volume filled at the respective time step. The reason for using the fill time sensitivity 
with respect to viscosity is that the fill time is a direct function of viscosity rather than 
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inlet temperature. This fact is reflected in the calculations performed, where viscosity is 
the only filling parameter calculated with the results from the thermal analysis. The new 
method described by equation (66) for estimating the fill times shows reduced error with 
the modification presented when compared to the method from equation (65), for multiple 
temperature ranges. The results using equation (65) and equation (66) are shown graphically 
in Figures 7, 8, and 9 for inlet temperatures of 300°, 350°, and 420°, respectively. In the 
legends, T in i et refers to the Taylor series initially used for the error estimation, and Visco. 
and T in i e t refers to the modified Taylor series estimation. The lines representing the viscosity 
based estimations more closely follows the finite element results than the lines representing 
the Taylor series estimation with inlet temperature only. It is observed that the viscosity and 
inlet temperature line does not always fall below the finite element results, i.e., tangent to 
the results curve. This error is explained by the fact that the filling calculation is dependent 
on the resin viscosity, which is in turn dependent on the calculated temperature profile. 



Figure 7: Comparison of the two methods for estimating fill time for T in i et = 300.0°F, for 
the axi-symmetric verification model. 
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Figure 8: Comparison of the two methods for estimating fill time for T inlet — 350.0°F, for 


the axi-symmetric verification model. 



Figure 9: Comparison of the two methods for estimating fill time for T in i et = 420.0°F, for 
the axi-symmetric verification model. 

4. Concluding Remarks 


In this report, the continuous sensitivity equation has been developed for non-isothermal 
considerations in RTM process modeling. This includes the non-isothermal RTM filling CSE, 
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the temperature CSE, and the cure CSE. Sensitivity parameters have been examined for the 
cases of material properties, boundary conditions, and geometric parameters. Numerical 
results were presented for an example axi-symmetric model and discussed. Finally, the fill 
time sensitivity results for the non-isothermal resin transfer molding numerical model were 
verified. This was accomplished by modifying the Taylor series initially used to estimate the 
fill time for small changes in inlet temperature as given in equation (65), with the modified 
Taylor series given in equation (66). The results using the modified Taylor series indicate 
that the dependence of fill time on inlet temperature is not direct but rather is a side effect 
of changes in viscosity. 
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